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ABSTRACT 

SN1998bw is the most luminous radio supernova ever observed. Previous discussions 
argued that its exceptional radio luminosity, ~ 4 x 10^® erg s“^, must originate from 
a highly relativistic shock which is fully decoupled from the supernova ejecta. Here 
we present an alternative model in which the radio emission originates from a 
sub-relativistic shock with a velocity ~ 0.3c, generated in the surrounding gas by the 
expanding ejecta. In this model, thermal electrons heated by the shock to a relativistic 
temperature of ~ 60 MeV, emit synchrotron self-absorbed radiation in the post-shock 
magnetic field. This model provides an excellent fit to the observed spectra provided 
that the thermal electrons are in equipartition with the ions behind the shock. The 
required magnetic field is much weaker than its equipartition value and could have 
been carried out by the progenitor’s wind prior to the supernova explosion. According 
to this model, the radio emission from SN1998bw is unrelated to the highly relativistic 
blast wave that produced the 7 -ray burst GRB980425. 


Subject headings: supernovae: individual (SN1998bw) - radio continuum: general - 
gamma rays: bursts 


1. Introduction 

The optical emission spectrum of SN1998bw classifies this event as a Type Ib/Ic supernova 
(Lidman et al. 1998; Sadler et al. 1998; Patat &: Piemonte 1998; Galama et al. 1998), suggesting 
that it resulted from a core collapse of a massive star which lost its hydrogen and helium envelope. 
Radio monitoring of the source revealed an exceptionally bright emission which peaked after 
12 ± 2 days (Kulkarni et al. 1998). The inferred radio luminosity at the source redshift of 0.0083, 
~ 4 X 10^® erg s“^, is the highest ever observed for a supernova. Multi-frequency measurements 
around the peak of the lightcurve showed a strong suppression of the emitted flux below a 
frequency of a few GHz. 

Aside from being associated with a rare class of supernovae, SN1998bw also occurred inside 
the 8 arcminute error circle of the Gamma-Ray Burst GRB980425 which was detected by the 
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Beppo-SAX satellite at about the same time (Sofhtta et al. 1998). The small likelihood for a 
chance overlap between the two events, ~ 10“^, led to the suggestion that the two might be 
associated (Galama et al. 1998) and to the conjecture that perhaps all SN Ib,c events lead to 
Gamma-Ray Bursts (Wang &: Wheeler 1998; Woosley, Eastman, & Schmidt 1998). However, the 
non-Euclidean number-count statistics, the event rate, and the temporal and spectral properties 
of most GRBs imply that SN1998bw-like events could only be associated with a minority fraction, 
^ 10% of the GRB population (Bloom et al. 1998). Indeed, the lack of an X-ray afterglow at the 
position of SN1998bw makes its potential GRB counterpart rare, since all but one of the other 
11 GRBs detected so far by Beppo-SAX showed an X-ray afterglow. In fact, a separate X-ray 
transient source, 1SAXJ1935.3-5252, was identified inside the Beppo-SAX error circle, with a flux 
comparable to that of the X-ray afterglows of other GRBs. This source provides an alternative to 
the SN-GRB association. A chance overlap between the error circle of GRB980425 and SN1998bw 
would fit better theoretical expectations, since supernova models have difficulties accounting 
for the highly-relativistic shock required to produce a GRB (Woosley et al. 1998). Hence, in 
gauging the likelihood of an association between SN1998bw and GRB980425 it is of fundamental 
importance to understand whether the bright radio emission from SN1998bw requires a relativistic 
shock by itself. If it does, then a more relativistic incarnation of the same shock at earlier times 
would serve as a plausible source for GRB980425 and hence strengthen the case for an association 
between the two events. 

The velocity of the supernova ejecta was measured from the blue wing of the Ca H line to be 
~ 60,000 km s“^ = 0.2c about a week after the explosion (Stathakis et al. 1998; Kulkarni et al. 
1998). Kulkarni et al. (1998) argued that the radio emission must have originated from a shock 
which was fully decoupled from the main supernova ejecta and expanded at a relativistic speed. 
They based their assertion on two arguments; 

1. If the radio source expanded at the measured ejecta velocity, then the brightness temperature 
of the source exceeded the threshold temperature for the so-called “inverse Gompton catastrophe” 
(Kellerman & Pauliny-Toth 1969, Readhead 1994). Under these conditions, Compton scattering 
of the radio photons to higher energies would have dominated the luminosity and violated upper 
limits on the X-ray flux of SN1998bw (Plan et al. 1998a,b). Relativistic expansion alleviates this 
problem. 

2. The lack of a strong variability at low frequencies implies that the source size exceeded the 
refractive scintillation scale of ~ 10^® cm at the peak of the lightcurve, and therefore expanded at 
a speed 0.3c. 

In order to put these constraints in a concrete physical context, we have constructed the 
simplest model for synchrotron emission due to the sub-relativistic ejecta of SN1998bw. In this 
model, the ejecta push as a piston and generate a strong forward shock in the surrounding gas 
(the latter being a remnant of the progenitor’s wind prior to the supernova explosion). The shock 
front heats electrons to a relativistic temperature, and these emit incoherent synchrotron radiation 
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in the post-shock magnetic field. To our surprise, we had found this model to be consistent with 
all observational data. In particular, the two constraints mentioned above are satisfied for the 
following reasons; 

1. The relativistic electrons responsible for the observed synchrotron radiation have a typical 
Lorentz factor of only ~ 360, and hence scatter the radio photons up in energy only into the 
optical regime but not into the X-ray band. Hence, the upper limit set by Beppo-SAX on the 
X-ray luminosity of SN1998bw is not in conflict with this model. In fact, the X-ray luminosity is 
dominated by scattering of the optical supernova photons by the relativistic electrons. This flux 
is consistent with the BeppoSAX bound for an expansion speed ^ 0.3c or a corresponding bulk 
Lorentz factor of only ^ 1.05. 

2. We find that a shock velocity of ~ 0.3c is sufficient for explaining the radio properties of 

SN1998bw. This velocity is ~ 50% higher than the measured velocity shift of the Call line. 
However, since the energy (~ 10^® erg) and mass (~ associated with the radio emitting 

plasma are small fractions of the total energy and mass of the ejecta, the required shock can be 
produced by a small amount of material preceding the bulk of the ejecta. In fact, our results are 
consistent with the modeling of the optical emission from SN1998bw by Woosley et al. (1998, see 
their §3.2), who suggested that the outer Mq of the supernova ejecta might have moved at a 
speed of ~ c/3 and carried r\j IqSO ergs. 

In section §2 we describe our model and derive the expected emission spectrum. In §3 we 
compare model predictions with the radio data, and derive the physical properties of the emitting 
electrons. The implications of our results are discussed in §4. Throughout the discussion, we adpot 
a distance of d = 38 Mpc to SN1998bw, assuming that it is located in the galaxy ESO 184-G82 at 
a redshift z = 0.0083 (Lidman et al. 1998; Tinney et al. 1998) and that the Hubble constant is 65 
km s“^ Mpc“^. 


2. Model 

In our model, the radio luminosity of SN1998bw is emitted by the material behind a forward 
shock, which is generated by the expanding ejecta in the medium surrounding the supernova 
progenitor. Since the ejecta expand at a sub-relativistic velocity, we assume to leading order that 
the radio flux observed at any given time is emitted by a static spherical shell of radius r. The 
compressed material behind the shock is expected to occupy a thin shell of width r/r], with rj ~ 10. 
We denote the mean number density of electrons inside this shell by n. 

The energy distribution of radiating electrons can be constrained from the observed 
synchrotron spectrum of the source. Kulkarni et al. (1998) have assumed that the radiating 
electrons have a power-law energy distribution. The synchrotron self-absorption spectrum of such 
an electron population is expected to show a spectral flux index of 2.5 at low photon frequencies 
(Rybicki & Lightman 1979). However, as noted by Kulkarni et al. (1998), the observed spectral 
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index is closer to ~ 2. This indicates that the electron distribution is truncated at low energies, 
as would be the case if they were thermal. For simplicity, we assume a mono-energetic electron 
population which should mimic closely the emission properties of thermal electrons. As shown 
below (see Figs. 1 and 2), our simple model provides an excellent fit to the observed spectra. 


The synchrotron power per unit frequency v emitted by a single electron of Lorentz factor 7 
is given by (Rybicki & Lightman 1979), 


P{v,^) = 
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where e and m-e are the electron charge and mass, B is the post-shock magnetic field strength and 
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The function F{x) describes the synchrotron power spectrum (Rybicki & Lightman 1979), 
averaged over an isotropic distribution of pitch angles. The total number of electrons in the 
radiating shell is n x (dvrr^/r/), and so the total flux per unit frequency observed at a distance 
d = 38 Mpc is. 
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where x (r/r]) is the optical depth per unit frequency for synchrotron self-absorption 

across the shell thickness. The self-absorption coefficient is given by (Rybicki & Lightman 1979), 
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For mono-energetic electrons of Lorentz factor 7 , dn/d^' = n5{'^' — 7 ), and the absorption 
coefficient can be integrated by parts. We then get, 
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where the functions, 
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The shape of the model spectrum, i{iy/i'c), is determined by a single dimensionless parameter, Tc- 
The frequency and flux normalizations are determined by the dimensional parameters t'c and A. 
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3. Comparison with Observations 

The most stringent constraints on the model parameters are obtained, as shown below, from 
observations near the peak of the radio lightcurve, ~ 12 days after GRB980425. These constraints 
are derived in §3.1. Observations at later times are discussed in §3.2. 


3.1. Observations Near the Peak of the Radio Lightcurve 

Figure 1 shows the observational data from the first two epochs (12 and 15 days after 
GRB980425) when multi-frequency measurements were taken around the peak of the radio 
lightcurve (Kulkarni et al. 1998). Our simple model of emission from mono-energetic electrons 
(solid line) provides an excellent fit to the data. The spectral shapes at both epochs are consistent 
with the model spectral shape in equation (P) for Tc = 0.6. The frequency and flux normalizations 
are Uc = 7-6 GHz, A = 3.3 x 10“^^ g at the first epoch, and = 6.2 GHz, A = 4.0 x 10“^^ g at 
the second epoch. For comparison, we also show the predicted spectra for a power-law population 
of electrons, dn/d'y oc 7 “^, with p = 2 (dashed line) or p = 3 (dotted), which are not consistent 
with the data. Note that the shape of the low frequency tail is independent of all parameters in 
both models. 

The free variables in our model are the density and temperature (Lorentz factor) of 
the radiating electrons, and the magnetic field strength behind the shock. The spectral 
constraints 12 days after GRB980425 can be parametrized in terms of the radius of the emitting 
shell, ri 6 = (r/10^® cm), the synchrotron frequency, t'g = (^c/ 8 ) GHz, the flux amplitude, 

H 3 = (H/3 X 10“'^^ g), and the optical depth coefficient r.e = (rc/0.6). We adopt a value of 7 = 10, 
and note that rig = 1 corresponds to an average expansion speed of ~ 0.3c. 

Equation (^) yields the electron Lorentz factor, 

7 PS 3.6 X 10^A^r^i, (10) 

which is equivalent to an electron temperature of kT = ~ 60A^t^q MeV. The magnetic 

field strength can then be obtained from equation (^), 

H PS 2.2 X 10“2i^8^3 G. (11) 

Substitution of these results into equation (|9|) yields, 

n PS 5.3 X cm“^. (12) 

The inferred number density of radiating electrons is much higher than that of the interstellar 
medium and is likely to be associated with a wind that originated from the supernova progenitor 
before or during the explosion. The existence of a dense ambient medium is natural given that 
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the progenitor might have lost its hydrogen and helium envelope prior to the explosion (Patat & 
Piemonte 1998). 

The mass of ions associated with the radiating electrons is small, 


Mion ^ 1.1 X 10 ^Mq, 


(13) 


where ^ = 2^X2 is the atomic mass per free electron. For fully-ionized metal-rich material, 1 x 2 ~ 1- 
The total shock energy associated with this mass at a post-shock fluid speed V = O.SFs x c is 
given by 

£^shock ~ MionV^ = 1.7 X 10"^^fi2VlusT6AlriQ ergs. (14) 

This amounts to less than a percent of the total hydrodynamic energy carried by the supernova 
ejecta. The radiating shell contains a total electron energy of 
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Our model predicts that the electrons and ions are nearly in equipartition, 
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This is consistent with our model assumption that the radio emission is produced by thermal 
electrons. The post-shock magnetic field is several orders of magnitude weaker than its 
equipartition amplitude. 


Let us now consider the constraints imposed by the X-ray observations of BeppoSAX. The 
radio radiation energy inside the shell volume is 

L^radio ~ = 5,3 X lO^^rie ergs, (18) 

where the factor of 4 is due to the fact that synchrotron radiation is approximately isotropic in the 
shell. Thus, the ratio between the synchrotron luminosity and the inverse Compton luminosity 
due to scattering of radio photons is 
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(19) 


A single Compton scattering boosts the synchrotron photon frequency up to a value of 
~ (4/3)7^t'c ~ 1-2 X Hz, i.e. into the optical regime. The resulting optical luminosity 

due to the upscattered synchrotron radiation is well below the thermal supernova emission in this 
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band - which is ~ 4 orders of magnitude larger than the radio luminosity. The upper limit on the 
X-ray flux of SN1998bw in the 1-10 keV band (Pian et ah 1998a,b) is obviously in accordance 
with the negligible X-ray flux due to the scattering of radio photons. 

Inverse-Compton scattering of optical supernova photons by the relativistic electrons yields 
a 7 -ray flux which peaks at a photon energy of ~ 300 keV and has a low-energy tail extending 
into the soft X-ray range. For input radiation of flux Fq at a frequency i'q, the low-energy tail 
of the inverse-Compton emission due to electrons with a Lorentz factor 7 is given by (Rybicki & 
Lightman 1979), fi, = where tt is the Thomson optical depth of the shell. 

For the supernova optical emission we approximate Fq/uq by the peak flux of ~ 10 mJy in the 
V-band {vq = 5 x lO^^Hz) and obtain 

Fx erg cm"^ s“^ ( 20 ) 


for the 1-lOkeV flux. The predicted flux is consistent with the BeppoSAX upper limit, 

~ 10“^^erg cm“^ s“^, for rig ^ 1. (Note that the BeppoSAX limit applies to day 8 , when both 
the optical and the radio fluxes are smaller by a factor ~ 2 compared to the corresponding 
fluxes at day 12 considered here). Inverse-Compton scattering of far infrared photons with initial 
wavelength ~ 30/um, produces a flux that peaks in the X-ray band. However, since the supernova 
flux at low frequencies is roughly thermal, fy (x the low-frequency tail of the inverse-Compton 
scattering of optical photons dominates the X-ray flux. 


Finally, we calculate the electron cooling time. The synchrotron cooling time is given by, 

bvrmpC 
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( 21 ) 


where cjt is the Thomson cross-section. The inverse-Compton cooling time equals, 

Fb 
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where, Fopt = Fopt^'/hc, is the optical radiation energy inside the shell volume, and we have used 
Fopt = 10^^erg s“^. Inverse-Compton cooling might therefore be responsible for the decline of the 
radio lightcurve immediately following the first peak. 


3.2. Observations After the Peak of the Radio Lightcurve 

Figure 2 shows the observational data around the second, weaker, radio peak (33 days after 
GRB980425) and at the last observational epoch (58 days after GRB980425). At both epochs, 
the spectrum is well htted by synchrotron emission of mono-energetic electrons, with a very small 
self-absorption optical depth, Tc = 0.06 and Tc = 0 at the earlier and later times, respectively. 
Since Tc <C 1, the spectrum is essentially fitted by an unabsorbed synchrotron model and for a 
given r the observations provide only two constraints on the three model parameters n, 7 and 



B. Thus, the late time observations can not determine the model parameters in the same way as 
observations near the first peak do (see §3.1). The decline in the self-absorption optical depth with 
time may be due to a decrease in the ambient gas density with radius, as expected for a wind. 
Similarly, the complicated lightcurve structure might reflect dumpiness in the spatial distribution 
of the ambient gas. 


4. Discussion 

The emission spectrum predicted by our model for mono-energetic electrons provides an 
excellent fit to the observational radio data of SN1998bw (see Figs. 1 and 2). This agreement 
implies that most of the radiating electrons share the same energy, as expected if their distribution 
is thermal rather than power-law in energy. Comparison of model predictions to data around 
the peak of the radio lightcurve suggests that the total energy content of the radiating electrons 
is comparable to the kinetic energy carried by the ions in the forward shock, implying that the 
two species are in equipartition at a temperature of ~ 60 MeV. On the other hand, the inferred 
magnetic field is much weaker than its equipartition amplitude and could have originated from a 
frozen-in seed field that evolved adiabatically with the ambient gas and did not suffer any dynamo 
amplification behind the shock. The field is likely to have been carried-out with the wind that 
emerged from the progenitor’s envelope prior to the supernova. In fact, its predicted amplitude of 
~ 20mG [Eq. ([TTD] is only an order of magnitude higher than that expected from the adiabatic 
compression of the interstellar ~ fj,G field as the interstellar medium density is increased by ~ 5 
orders of magnitude to the value of the ambient gas density [Eq. (0)i- 

The source radius required in order that our model be consistent with X-ray observations 
near the peak of the radio emission, ~ 12 days after GRB980425, is only ~ 10^® cm. This implies 
a shock velocity of ~ 0.3c and a corresponding Lorentz factor of only 1.05. For these parameters, 
the mass and energy associated with the radio emitting plasma are ~ 10~^Mq and ~ 10^® ergs, 
respectively. Our model is therefore consistent with the modeling of the optical emission by 
Woosley et al. (1998) who suggested that the outer 10~^Mq of the ejecta in SN1998bw might 
have moved at a speed of ~ c/3 and carried r\j 1050 gpgg^ 

Based on these findings we conclude that the radio luminosity of 1998bw could have been 
produced by the sub-relativistic shock in front of the supernova ejecta. In our model, there is 
no need to invoke an ultra-relativistic shock that might have been related to GRB980425. This 
raises the possibility that GRB980425 was associated instead with the transient X-ray source 
1SAXJ1935.3-5252. 

The large scatter in the radio luminosities of Type Ib/Ic supernovae (for examples, see Weiler 
&: Sramek 1988; van Dyk et al. 1993; and Weiler et al. 1998) is not surprising in view of the fact 
that the total energy and mass of the radio-emitting plasma are negligible fractions of the overall 
energy and mass budgets of the supernova ejecta. The scatter might also be enhanced by viewing 
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angle variations due to non-spherical explosion geometries (Hoflich, Wheeler, & Wang 1998). 

Our discussion ignored various effects, such as the increase in the self-absorption optical depth 
towards the limb of the source, the free-free absorption in the surrounding stellar wind, or the 
nonzero time delay due to the finite light crossing time across the source. Future work is necessary 
in order to refine the detailed predictions of this model. 
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Fig. 1.— Emission spectrum from SN1998bw around the peak of the radio lightcurve. The 
data points from Kulkarni et al. (1998) are shown 12 days (squares) and 15 days (circles) after 
GRB980425 and are compared to theoretical predictions for a population of electrons which is 
either mono-energetic (solid line), or power-law in energy with a spectral index of p = 2 (dashed) 
and p = 3 (dotted). The frequency and flux of the second epoch observations (15 days after 
GRB980425) are renormalized so as to show that the spectra at both epochs have a similar shape 
(frequency multiplied by 1.23, flux by 1.24). The model parameters for mono-energetic electrons 
are Tc = 0.56, A = 3.3 x 10“^^ g, Vc = 7.6 GHz [see Eq. (^]. 
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Fig. 2.— Emission spectrum from SN1998bw around the second radio peak (circles) and at the 
last observational epoch (squares), 33 and 58 days after GRB980425 respectively (Kulkarni et ah 
1998). The data points are compared to the theoretical predictions for a mono-energetic electron 
distribution with Tc = 0.06 at the first epoch and Tc = 0 at the second epoch. 




